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The paper is devoted to the introduction of simple analytical relationships between statistical distributions of 
various radiative transfer characteristics for an inhomogeneous turbid layer with the extinction coefficient 
varying in the horizontal direction. Results are valid for an optically thick light-scattering layer having arbi- 
trary local scattering laws and single-scattering albedos. It is shown that the statistical distribution of the 
optical thickness can be obtained directly from the measured statistical distribution of the reflectance or trans- 
mittance of a horizontally inhomogeneous light-scattering layer. © 2005 Optical Society of America 
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1. INTRODUCTION 


Radiative transfer characteristics of horizontally inhomo- 
geneous media (e.g., clouds and snow) can be studied us- 
ing Monte Carlo numerical calculations.” This technique 
allows one to obtain accurate results for arbitrary spatial 
distributions of scatterers. However, for such calculations 
to be exact, a lot of computational time for optically thick 
media (especially at wavelengths where light absorption 
is weak) is required. Also, simple physics behind relation- 
ships between physical parameters under study is hidden 
in numerical results. 

Kokhanovsky? established relationships between vari- 
ous radiative transfer characteristics of horizontally inho- 
mogeneous turbid media (e.g., the average diffuse trans- 
mittance and the average bidirectional reflection 
function) in the framework of the asymptotic theory as 
given in Refs. 4 and 5. The results obtained are valid for 
optically thick weakly absorbing media only. 

The task of this work is to extend the study of 
Kokhanovsky’ to the case of arbitrary absorption. An- 
other difference from the previous work is that direct re- 
lationships between statistical distributions of various op- 
tical characteristics are derived in the case of thick turbid 
layers. The previous work was concentrated on average 
characteristics and coefficients of variances. 

The limitation with respect to the optical thickness 7 of 
a turbid layer remains. So the technique cannot be ap- 
plied to the case of optically thin turbid layers. 


2. OPTICAL PROPERTIES OF 
INHOMOGENEOUS SCATTERING MEDIA 
OUTSIDE ABSORPTION BANDS 


We assume that scatterers are confined in a plane- 
parallel slab illuminated from above by an infinitely wide 
uniform light beam (e.g., clouds illuminated by the Sun). 
The angular distribution of light scattered by an elemen- 
tary volume of a turbid medium is assumed to be constant 
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at any given place inside the medium under study. Effects 
of light absorption are neglected, which is an accurate as- 
sumption in the spectral regions positioned far from ab- 
sorption bands of substances present in the scattering 
layer under study. However, we assume that the extinc- 
tion coefficient fluctuates in the horizontal direction (e.g., 
owing to fluctuations in the concentration of scatterers). 
Such fluctuations lead to fluctuations of the optical thick- 
ness. Strictly speaking, this makes one-dimensional (1-D) 
radiative transfer theory® (RTT) not applicable to the case 
under study. 

However, a slight modification of the theory allows us 
to extend the applicability of the 1-D RTT to treat (al- 
though approximately) this more complex case. In par- 
ticular, it is assumed that the reflection function and also 
other radiative characteristics of a slab can be presented 
as integrals: 


F- | Ynf(aar, (1) 


0 


where Y is a radiative transfer characteristic under study 
and f,(7) is the statistical distribution of the optical thick- 
ness 7. This approximation is called the independent pixel 
approximation. The accuracy of this approximation as 
compared with Monte Carlo calculations was studied in 
Refs. 2, 7, and 8. Barker’ proposed the use of the gamma 
distribution 


Fera) 
= ———— 7 exp| - u— 2 
f7) ATD ss a (2) 


in conjunction with the two-stream approximations’? for 
the characteristics Y(7) to account for the horizontal inho- 
mogeneity effects of turbid layers. Then the integral in 
Eq. (1) can be evaluated analytically. The statistical dis- 
tribution given by Eq. (2) is shown in Fig. 1 at u=8 (the 
coefficient of variance equals 1/3) and 7)=20. 
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Fig. 1. Statistical distribution of the optical thickness calcu- 
lated using Eq. (2) at 7)=20 and u=8. 


Such a technique leads to simple analytical equations 
for a number of radiative transfer characteristics impor- 
tant in global circulation models. However, the two- 
stream approximation cannot be used for the bidirectional 
reflectance and transmittance calculations. Then the 
asymptotic RTT can be used.’ It follows in the framework 
of this theory”! for a homogeneous optically thick non- 
absorbing plane-parallel turbid layer that 


R¢é, N, ¢) = Ro..(€, N, a) = tKo(€)Ko(), (3) 


T(é, n) = tKo(2)Ko(), (4) 


where R(é,7,¢) is the reflection function, T(é, 7) is the 
transmission function, Ro..(é,7,¢) is the reflection func- 
tion of a semi-infinite nonabsorbing medium having the 
same optical properties as a finite layer under study, Ko(é) 
is the escape function for a nonabsorbing layer deter- 
mined via Ro..(é, 7, ¢) ast 


3 1 20 1 
Kl) =7 etf af dole) t. (5) 
TJo 0 


é is the cosine of the zenith illumination angle and 7 is 
the cosine of the zenith observation angle, ¢ is the rela- 
tive azimuth, and 


1 


t= 
a+br 


(6) 


is the global transmittance defined as 
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and we accounted for the fact that the light diffused 
through thick layers is azimuthally independent [see Eq. 
(4)]. The parameters a and b are given as follows”: a 
=1.07 and b=0.75(1-g) for arbitrary phase functions. 
Here g is the asymmetry parameter. This parameter is 
derived from Mie theory in the case of spherical 
scatterers.!* For instance, its value is close to 0.85 for wa- 
ter clouds in the visible.” 
The escape function is normalized as follows": 


1 
af Ko() nd y= 1. (8) 


0 


Sobolev" proposed the following approximation for the 
escape function valid at arbitrary phase functions and 7 
>0.2:Ko(7)=3(1+27)/7. This relationship allows one to 
derive simple equations for the diffuse transmittance, 


1 
wo=2f T(é,n) ndn, (9) 
0 
plane albedo, 
1 
rao =2 R(é,) nn, (10) 
0 
and the spherical albedo, 
1 
r=2{ raléédé, (11) 
0 
where 
B 1 2m 
R(é, n) = >Í R(é, 7, ¢)de (12) 
27 J 9 


and similar for T. So we have!® 


ta) =tKo(g), ralé) = 1-tK) (8), 


where we accounted for the fact that rg..=1, r.=1 by defi- 
nition for nonabsorbing semi-infinite media. 

Radiative characteristics r, t, ralé), tg(6), R(é, 7, ¢), and 
T(é, 7) can be measured. In particular, pyranometers in- 
stalled on an aircraft are used to measure r,(é) over a 
cloud field. 

We pose the following question. Is it possible to obtain 
the statistical distribution of the optical thickness /,(7) 
from the measurements of the statistical distribution of 
any radiative transfer characteristic fy(Y) as specified 
above (Y=r,t,...), provided that conditions valid for the 
theory presented here are fulfilled? The answer on this 
question is positive. 

To show this, we use the following relationship derived 
in the probability theory”: 


r=1-t, (13) 


fy(Y) = f7). (14) 


7 
dY 
This formula is valid for monotonic (increasing or de- 
creasing) functions Y(7). Radiative characteristics in the 
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asymptotic regime (7+) satisfy this condition." There- 
fore, we can use Eq. (14) to relate fy(Y) to f(7) in conjunc- 
tion with Eqs. (3), (4), (6), and (13). In particular, we ob- 
tain from Eqs. (6) and (14) (see Fig. 2) 


Ft) =DE FAT, (15) 


where r=b7!(t-!-a) [see Eq. (6)]. This means that f,(7) 
can be obtained from the measured function f;(t) as fol- 
lows: 


FAD = bt (af ,(E(0). (16) 


Functions T(é, 7) and t,(é) are related to ¢ linearly. This 
means that 


f(T) =f)Ko' AK), — fi (ta) =f) Ko'(8), 


f(T) =f ta) Kon), (17) 
and, therefore, 


FAD) = DPF DKE, FAD) = bt? f(T) Ko(€) Koln). 


(18) 


We also find that functions fr(R), f- Ta)» and f,(r) coincide 
with f(T) ,f;,(ra),f:(), respectively [see Eqs. (3) and (13)]. 
Therefore, one concludes that measurements of the statis- 
tical distributions of radiative characterististics above 
and below a scattering nonabsorbing turbid layer can be 
used to derive the optical thickness distribution using the 
analytical equations given above. The results are quite 
general and not restricted to specific types of statistical 
distributions. 
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Fig. 2. Statistical distribution of the total light transmittance 
by a turbid layer having the statistical distribution of 7 shown in 
Fig. 1. Calculations were performed using Eq. (15) at g=0.85 and 
assuming that r=(t-!—a)b~! [see Eq. (6)], where b=0.75(1-g). 
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3. ACCOUNTING FOR LIGHT ABSORPTION 


All substances at all wavelengths absorb some portion of 
incident light for real-world scattering media. Therefore, 
it is of importance to study how the effect of light absorp- 
tion in a turbid layer influences the results given above. 

The starting point is the generalized form of Eqs. (3) 
and (4)” for a homogeneous light-scattering plane- 
parallel layer valid as 7— œ for arbitrary absorption and 
phase function: 


R(Eé, 7,9) =RAE, 7,9) - TE Ne, (19) 
T(E, n) =tC°K(K(n), (20) 
where 
MC? e kt 
a ra = 


is the global transmittance for an absorbing turbid layer. 
Functions R..(é,7,¢) and K(ņ) have the same sense as 
Ro.(é,7,¢) and Ko(7), respectively, but for an absorbing 
layer. The constant C is defined as 


1 
c=2f Kn) ndn, (22) 


0 


which is equal to 1 for nonabsorbing layers [see Eq. (8)]. 
The procedures to calculate the auxiliary functions 
R.(é,7,¢) and K(ņ) and constants M,N,k are well 
known.!""}3 It is of importance for us that these functions 
and constants do not depend on the optical thickness. 
They are determined solely by the phase function p(@) and 
the single-scattering albedo wo. These characteristics are 
assumed to be constant inside a turbid layer under study. 
In particular, it follows as wọ—> 1 (Ref. 13): k ~[38(1- wọ) 


x(1-8)]?,  MK(K(»)~(1-exp(-2y))Ko(Kol), Rx 
=Ro.. exp(-yRo1Ko(é) Ko(), N~exp(-y), where y 
=4,(1-w)/3(1-g). Useful approximations for 


Ro.(é,7,¢) can be found in Ref. 13. 

Equations (19)-(21) give the dependence of reflection 
and transmission functions with respect to 7 in the ana- 
lytical form. So they can be used to find the derivative as 
specified in Eq. (14) analytically. As in the nonabsorbing 
case, we start from the search of the relationship between 
f(t) and f,(7). 

For this, we express 7 via t using Eq. (21). This gives 


| 2tN? | 
tr=k Inj —~——_—-/,, (23) 
MC*(B(t) - 1) 


AN??? 
Bit) = 1+ wc!’ (24) 


Therefore, we have for the derivative 


dr 1 4N??? 
== i= : (25) 
dt kt B(B -1)M?c+ 


where 


where ¢ is given by Eq. (21), and we omitted arguments. 
This equation is simplified as £— 0: 
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dr 1 
an (26) 
Therefore, introducing 
1 AN??? 
oe al "- ME- “ah i 
we obtain 
fit) = O(n f(7). (28) 


This analytical equation allows one to determine the sta- 
tistical distribution of the optical thickness from the glo- 
bal transmission measurements [see Eq. (7)] as 


f(a) = PDRE). (29) 


Using Eqs. (9), (20), and (22), we derive for the diffuse 
transmittance 


ta(é) =tC'K(é). (30) 
This means that 
fita) = CK HEF). (31) 
Similarly, one finds 
fr(T) = CKK fA). (32) 


So we have 


FAD = DACKO, (ta); 


FAD = PAC PKOK fT). (33) 


It follows that the statistical distributions of the light 
transmittance under different measurement schemes 
(e.g., diffused or direct light illumination conditions) are 
interrelated and governed by the statistical distribution 
of the optical thickness. Equations (29) and (33) can be 
used to determine the statistical distribution of the opti- 
cal thickness from measurements of the transmitted 
light. 

Let us consider the reflectance now. First of all, we 
write, taking into account Eqs. (10), (11), and (19), 


ral = Taal) - tCOK(EN e7 (34) 
and 
r =r- tN e™%7 (35) 
or 
MNC? eo 2k 
rafe= Tonto (36) 
So we have from Eq. (36) 
1 MNC? 
r= — In) N? + 3 (37) 
2k ra -r 


dr 
rn =[2R(r..-—r)(1+M4NC~(r.. —1r))}1. (38) 
r 


Therefore, introducing 
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Q(r) =[2k(r2.-r)1+MINC~(r.,-7)) Tt, (39) 
we have 


fAr)=Q)FA7D, FAD = QIO), (40) 
where r(7) is given by Eq. (36). This equation allows us to 
find the statistical distribution of the optical thickness 
from measurements of the spherical albedo r. Similar re- 
sults can be obtained for fR(R),f- fa). Namely, we have, 
following the procedure given above, 


FAD =DR)fRE),  FAD=Wr Df ra), (41) 


where 


D(R) =[2k(r..-r)(1+M"NCPK(K(n\(r.-7r))T", 
(42) 


W(rq) =[2R(r..-r)(1 + M NCK Alr- r). (43) 


Equations (29), (33), (40), and (41) can be used to relate 
statistical distributions for reflectance and transmittance. 
Results for the statistical characteristics of light absorp- 
tion processes such as a=1-r-t and a,(é)=1-rg( 
—tg(é) can be also easily derived from the equations given 
above. 


4. CONCLUSION 


In conclusion, we underline that, as expected, reflection, 
transmission, and absorption statistical distributions are 
determined by the statistical distribution of the optical 
thickness of a random medium under study. All these dis- 
tributions are interrelated as clearly demonstrated in this 
work using the analytical approach valid as 7—~ for ar- 
bitrary absorption and phase functions. 

The statistical distributions derived can be used to find 
the average values of correspondent characteristics, their 
higher-order moments (e.g., standard deviations), and re- 
lationships between various moments (see, e.g., Ref. 16). 

The derivation of statistical distributions of the optical 
thickness and also correspondent distributions for trans- 
mittance from reflectance measurements are of impor- 
tance in a range of remote-sensing techniques (e.g., in 
cloud airborne and satellite remote sensing!™!8). 

The derived equations can be used also for the determi- 
nation of the optical thickness of a homogeneous random 
layer from reflection or transmission measurements [see, 
e.g., Eq. (37)]. However, the auxiliary functions and pa- 
rameters must be known in advance for such a procedure. 
Simple approximate equations for them [dependent on 
the pair (w9,g)] are given in Refs. 4, 5, 13, and 19. 
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